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Abstract: In this presentation, parallel multilevel characteristic basis function method 
(MLCBEM) has been hybridized with the randomized pseudo-skeleton approximation method 
(RSPA) and randomized singular value decomposition (tS VD) method, for the analysis of the 
scattering from electrically large rough surfaces at low grazing angles. MLCBFM defines the 
Characteristic Basis Functions (CBFs) on a larger domain and thereby achieves a higher 
compression rate. The reduced matrix has a much smaller size compared to that of the original 
system, which enables us to use a direct solver, rather than an iterative one. The RSPA algorithm 
accelerates the generation of the reduced matrix while the rSVD expedites the generation of 
CBFs. The hybrid method is found to be both accurate and efficient, and we use it in this work 
to investigate the problem of bistatic scattering from Gaussian and exponential rough surfaces 
at low grazing angles. 
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MOTIVATION (1) 


“Higher compression rate can be achieved by 
using larger domain in CBFM!1), However, the 
increased degree of freedom (DOF) in large 
domain makes the generation of the 
Characteristic Basis Functions (CBFs) very 
time- and memory-consuming. The reasons 
are: 1) Impedance matrix pertaining to each 
domain should be calculated; 2) Singular 
Value Decomposition (SVD) is adopted to 
remove the redundancy between the CBFs. 


“To mitigate such problem, Multilevel CBFM 
(MILCBFM) has been proposed. 


[R1] J. Laviada, F. Las-Heras, M. R. Pino and R. 
Mittra, "Solution of Electrically Large Problems With 
Multilevel Characteristic Basis Functions,” Antennas 
and Propagation, IEEE Transactions on, vol. 57, no. 
10, pp. 3189-3198, October 2009 
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MOTIVATION (2) 


“*Similar to MLCBFM, a hybrid approach, namely 
the hybrid CBFM/ACA/UV method, was 
developed in our recent work to enable the usage 
of large-size blocks [R2]. 


“+ However, the increased number of unknowns 
in large-size blocks leads to an increase in the 
computational cost, because it becomes 
necessary to carry out the Singular Value 
Decomposition (SVD) to remove the redundancy 
between the CBFs and, as a consequence, the 
overhead for generating the CBFs is dominated by 
the computational cost of the SVD operation. 


Ground 


Incident Field 


Random scattering 
[R2] Chao, L. and R. Mittra, “Accelerating solution of rough 4 
surface scattering problems by using the UV technique in 
conjunction with the Characteristic Basis Function Method 
and the Adaptive Cross Approximation ” IEEE International 
Symposium on Antennas and Propagation and USNC-URSI 


Radio Science Meeting, 2015. 
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CHARACTERISTIC BASIS FUNCTION METHOD 


Step - I: The first step in CBFM is to divide the 
problem in smaller regions with respect to the 
Original geometry, called “blocks” [R,] 


% The induced currents in each blocks are 
represented by using a type of high-level basis 
function defined Primary Basis Functions 





Step - Il: The Primary Basis Functions are ere 
evaluated by analyzing each block in the 
absence of the others: 


j =1,2,--,M 
7 +P, =V. 
=i Th —hnh y= 1,2,---,N 





[R3] E. Lucente, A. Monorchio and R. Mittra, “Parameterization of the tapered incident wave for numerical 
simulation of electromagnetic scattering from rough surface," Antennas and Propagation, IEEE Transactions on, 
vol. 56, no. 4, pp. 999-1006, April 2008 
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CHARACTERISTIC BASIS FUNCTION METHOD 


% Each block is excited by using a Spectrum Of Plane Wave Spectrum on Block 
Plane Waves (SPWs) intentionally overestimating : | 
the Degrees Of Freedoms (DOFs) of the problem 


Step - Ill: After generating the CBFs, a 
thresholding procedure based on the a Singular 
Value Decomposition (SVD) approach, can be 
employed in order to reduce the number of CBFs 
and the reduced matrix size 


Step - IV: The final current distribution J induced 
on the scatterer can be expressed as a linear 
combination of the CBFs provided by the SVD 
procedure (Ja m): 





n=1 n=l n=1 — 
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CHARACTERISTIC BASIS FUNCTION METHOD 


Step — V: The final step to be performed is the 
generation of the Reduced matrix Z" that can be 
accomplished by applying the Galerkin testing 
procedure employing the CBFs as testing 
functions 

J VJV 


i aA AEE 


IN 


Z: MoM impedance matrix 

J: matrix comprising all the CBFs after the SVD 
procedure 

V: RHS vector related to the real problem 
excitations 


*% The CBFM leads to a reduced matrix, which is 
much smaller than the original one, and this 
obviates the need for an iterative solution of 
problems requiring a large number of unknowns 
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Hybrid CBFM/UV/ACA Method 


%a hybrid approach, namely the hybrid CBEFM/ACA/UV method, was developed in our 
recent work to address the same issue. By using UV technique, the computational time 
and memory complexity required by the matrix filling process are decreased from 
O(B Nic )to O(B Neve log ac ), where B is the number of blocks and Wu. is the average 
number of RWG in each block. Comparing to MLCBFM, hybrid CBFM/ACA/UV has the 
advantage that it does not need to generate either the CBFs or the reduced matrix in 
the lower level (two-level domain decomposition policy is normally adopted for 
MLCBFM). 
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UV decomposition method: for scalar problem 


For scalar wave scattering, the rows and columns can be sampled uniformly in the 
source and field regions to select the most important interpolation points, and the 
method is called coarse-coarse sampling and proved to be quite efficiently . 


Region] Region2 


— 
= ae a = 
a 





Two well-separated regions with coarse-coarse sampling 


However, for vector wave scattering problem, where the RWG vectors basis function 
is employed, the EM interaction is not only related to the space distance but also 
affected by the current direction of the RWG function. Consequently, the matrix 
element varies abruptly and the smooth characteristic of the matrix elements that 
the original UV method demands doesn’t exist . 


*L.Tsang, Q. Li, and P. Xu, “Wave scattering with the UV multilevel partitioning method: 2. Three-dimensional 
problem of nonpenetrable surface scattering,” Radio Science, vol. 39, RS5011, 2004. 


New sampling method for vector problem 





0.1 
0 


50 100 150 200 
The EM interaction for a far interaction submatrix from the 3-D cube scattering problem. 
The distribution is disordered and oscillatory according to the natural space ordering. 





l 10 20 30 40 50 60 70 8 90 m 110 120 130 140 150 160 170 180 190 200 210 
ndex 

Re-ordeing of the EM interaction. It is observed, after sorting, the norm varies 

monotonically according to the new ordering, i.e. the interaction between the source 


and field regions are re-ordered by their strength from the strongest to weakest, form 
the point view of the physics. 


Acceleration of the matrix vector multiplication 


Matrix Vector Multiplication 


== N; — u = 
S == D (Zr) +) (Z jar), 
i=l j= 


= sub 


(Z i = ONT r<< min(m, n) 


) w= ~(n,) = ~(n,) 
=() (Znear);) LE” +() (Zar), I” 
i=l j=l 


= sub S) = Sn 


(Z far Joa Ima a U mxr (V rxn Twa ) 


Memory requirement and Computation complexity 
mxn == rx(m+n) 


O(N’) ===» O(NlogN) 
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Hybrid CBFM/UV/ACA Method 


% The UV technique is used to accelerate the fill-in of the submatrices pertaining to 
blocks in level 1 





% The ACA method is used to accelerate the generation of the reduced matrix in level 1 


Comparing to MLCBFM, hybrid CBFM/ACA/UV has the advantage that it does 
not need to generate either the CBFs or the reduced matrix in the lower level 
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CODE VALIDATION 


“* The developed code has been validated by analyzing the canonical example of a 
Perfect Electric Conductor (PEC) sphere 


% The results provided by the Hybrid 
CBFM/ACA/UV have been compared with 
the conventional CBFM ones and with the 
Mie solution 


% The code accuracy have been tested by 
using the reported Relative Error definition 


*% Sphere radius: 3.0m 
% Frequency: 300MHz 
% Number of RWGs: 52,023 








¢¢ SVD threshold: 1e3 RCS —RCS 
MCBFM ,n Mie,n 
% ACA threshold: 1e Rel. Error = 100* yy RSs Be 
n Mie,n 
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CODE VALIDATION 


Bistatic RCS: Performance summary: 


===Mie 
—— CBFM 
—--CBFM/ACA/UV 


CBFM CBFM/ACA/UV 


Number of blocks 


ofa? (dBsm) 


Dimensions of the 
reduced matrix 


Solution time(s): 





60 80 100 120 140 160 180 
Bistatic Angle (Degrees) 


Relative Error(%) = (2.3327 


“* Very good agreement has been achieved between the proposed approaches 
e Sphere radius: 3.0m 

e Frequency: 300MHz 

e Number of RWGs: 52,023 

e SVD threshold: 1e-3 

e ACA threshold: 1e-3 
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CODE VALIDATION 


% We have also tested the numerical efficiency of 
the parallel algorithm by progressively increasing 
the number of CPUs employed to solve the sphere 
problem 


“* The obtained efficiency has been compared with 
respect to the ideal case where the solution time 


scales as 1/n 
CPU Time: 


—*— Ideal Case 
—-@ —- Parallel MLCBFM 


Normalized CPU Time 


8 10 
Number of CPUs 


% The parallelized CBFM algorithm scales quite well 


Parallel Implementation: 









Multilevel Decomposition and Data Preprocessing 


Generating CBFs for each Block independently 


Processor Processor Processor 
l 


Using MPI Allgather to broadcast the CBFs 
Processor Processor Processor 


l 


Generating the Reduced Matrix 


Processor Processor Processor 
i 


Solve using 
LU Decomposition 
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THE randomized Singular Value Decomposition 


«* The randomized Singular Value Decomposition (rSVD) has been implemented in 
order to efficiently generate the CBFs of the large blocks in MLCBFM 


The rSVD begins by randomly projecting the original matrix to a lower-dimension 
matrix, while the range of the original matrix is asymptotically kept intact. The 
projected matrix, which is much smaller, is then factorized by using a full-matrix 
decomposition such as the SVD, after which the resulting singular vectors are back- 
projected to the original space. The algorithm can be readily implemented by using 
the following steps: 


1. Generate a collection of / random vectors of length N (from a Gaussian 
distribution) and arrange them in columns to from a rectangular matrix R 
with a dimension of N XI]. 

2. Compute an M X/ sample matrix Y=AR. The column vectors form a basis for 
the range of A. 

3. From an M X I matrix Q whose columns form an orthonormal basis for the 
columns of matrix Y. Then A=QQ°A. 

4. Form the /XN “small” matrix B= Q“A. Then A=QB. 

5. Form the SVD of B (cheap since B is “small”), B= uxv’. 

6. Form U=QU . 
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THE randomized Singular Value Decomposition 


The accuracy of the rSVD is considered by the test case of scattering by a 4AX4A 
rectangular plate, illuminated by a spectrum of plane waves incident from 
O° <0s180" andO <@s360° , with N,=18 and N,=36. We begin by calculating 
the coefficients of pre-SVD CBFs (characteristic basis functions) by using the 
Method of Moments (MoM), and then arrange the CBFs in columns to form a 
matrix J-BFs, which has a dimension of 1568 X 1296. Next, traditional SVD and rSVD 
are performed on J‘®*S. The post-SVD CBFs are derived and denoted as U>V? and 
UY? respectively. Following this we define a difference matrix U“'" as 


diff SVD SVD 
Ue =U" -U' 











15 . { x 10 

_ Traditional SVD 3 

0.9 : —&— randomized SVD 2 0.8 
A O 
= 

a 0.8 2 0.6 
> wo 
3 3 

2 0.7 o 0.4 
= a 
a © 

0.6 a 0.2 
p 
7 

0.5 0 

0 30 60 90 120 150 
Index 
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THE randomized Singular Value Decomposition 


Four pre-SVD CBFs matrices are generated by enlarging the dimension of a 
rectangular plate to show the efficiency of the rSVD. The dimensions are 
1568 X 1296, 2720 X 1296, 4608 X 1296 and 8368 X 1296. The number of random 
vectors / is chosen to be 550. The CPU times of the SVD and rSVD are shown in 
Table 1, from which we can see that rSVD is much more efficient than SVD, 
especially when the dimensions of the object are large. 


0.74s 7.63S 58.625 488.23s 


0.15s 0.35 0.94s 3.46s 
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Numerical Result-Rough Surface Scattering 


Two-dimensional Perfect Electric Conductor (PEC) rough surfaces are analysed by 
using MLCBFM. Three surfaces whose sizes are 8A X 8A(case 1), 16A X 16A(case 2) 
and 32A *32A(case 3) are considered. The root mean square (RMS) height is 
chosen to be 0.168A and the correlation length is assumed to be 10 times the 
RMS height. 100 realizations are implemented to derive the statistical results. 


----RS 8A by 8A 
----RS 161 by 16 À 
— RS 321 by 32 à 


Bistatic Coefficients (VV) [dB] 





3B -60 -30 0 30 60 90 
Scattering angles [deg] 


The angular widths of coherent bistatic scattering coefficients become smaller 
with an increase of the size of the surface 
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Numerical Result-Rough Surface Scattering 


As discussed above, we desire to keep the block sizes large enough to gain a 
high compression rate. Two different block sizes are adopted and we also use 
SVD to obtain the results, instead of the rSVD, for the sake of comparison. The 
computation time for a single realization is shown in Table 2. It is observed that 
when the block size is chosen to be large, using the rSVD in the MLCBFM 
algorithm results in significant time-saving. For the rough surface with 
dimensions of 32A X 32A, the total time saved for the derivation of the statistical 
result was on the order of 81 hours. 


“a 


size SVD rSVD SVD rSVD SVD rSVD 


2885S 243s x x 
251s 192s 1905s 1581s 8460s 6976s 
x 1678s 942s 7264s 4345s 
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Numerical Result-Radiation Problem 


~x Length Width Height 


156m 21m 28m 





Ship with a monople antenna 


% The operating frequency is chosen to be 
100 MHz and it requires 74,0262, 
unknowns to discretize the target. 
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Numerical Result-Radiation Problem 


% 32 Processors 


¢¢ 32 Blocks on the first-level 
% 128 Blocks on the second-level 

1790 *% Solve Time: 40 Mins 
% Peak Memory: 60GB(approximately) 
% The total CPU time is reduced from 
78 minutes to 46 minutes after 
combining MLCBFM with rSVD. 





180 
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Conclusion 


The rSVD is introduced in this work--for the first time in the context of CEM-- 
and combined with the MLCBFM to solve electromagnetic scattering 
problems using a domain decomposition approach. The proposed algorithm 
helps to efficiently derive the CBFs as well as to generate the reduced matrix 
when using large blocks. Numerical results for a statistical study of the 
problem of scattering from electrically large random rough surfaces are 
presented to show that the computation time is significantly reduced, in 
comparison to that needed in the conventional MLCBFM or Hybrid 
CBFM/ACA/UV, when the proposed algorithm is employed to solve large 
problems. 
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